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We find expectation values of functions of time integrated two-level telegraph noise. Expectation 
values of this noise are evaluated under simple control pulses. Both the Gaussian limit and 1/f noise 
are considered. We apply the results to a specific superconducting quantum computing example, 
which illustrates the use of this technique for calculating error probabilities. 



I. INTRODUCTION 

Two-level telegraph noise, sometimes called popcorn noise or burst noise, appears in a variety of sources. At any 
given time t, the random telegraph signal Y{t), which represents the derivative of the noise, is in either the positive 
state A or the negative state —A. It has probability ^ of flipping from the positive state to the negative state, and 

probabihty ^ of switching from the negative state to the positive state. After time t, the parameter of the noise is 

given by 6* = Y{t)dt. If there are no flips, then after time t the parameter will be ±0c, where 9c = At. An example 
of a Y{t) and its corresponding as a function of time is given in Fig. [TJ 




FIG. 1: An example Y{t) and 6 



In this paper, we show how to find the expectation value of a function f{0), which we shall write as E[f{9)]. We 
can write a function in the Fourier basis, and then find the expected values of E[e^"^^], which we show in Sec. |TT1 for 
the case in which the correlation times tq and ti are equal, and in the general case of tq 7^ ti in Sec. |Vl 

In Sec. mil we look at how control can be used to suppress effects of this noise on a pair of qubits. In Sec. IIVI the 
results are applied to a superconducting qubit system which looks promising for Quantum computing 



II. EVENLY DISTRIBUTED RANDOM TELEGRAPH NOISE 



Here we assume that the correlation time Tc = tq = ti is the same in both directions of flips. Suppose we want to 
find the expectation value E[f{6)], where 9 has a random telegraph distribution. 

Lemma II. 1 Since it is Poisson in the number of flips f, the distribution of 9 after time t is given by 

°° \/ 

=e-^^rf(0)/ — , for - 9c<9<9c, (1) 

where the distribution for a given number of flips is on this same domain 

( ±9c forf^O 
dfiO) = ^4;,,. ^ll]^ for odd f (2) 
{df_i{9) for even f > 
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Proof: First we look at the discrete case with q steps. g{h, c) = represents the number of ways to put c 

identical objects into b different boxes. We have r intervals where we are heading 1 to the right per unit of time, and 
I intervals where we are heading 1 to the left per unit of time. The probability of being at position j afterwards is 
the coefficient of in 

k)g{l,q-k) 

Since g{b, c) is asymptotically equivalent in c to j^rrpj, for large k and q, this becomes 



(r + /-l)! J^k'^-^q-k) 



l)!a- 1)! ^ 



/-I 



r+l-l 



-X 



2k-q 



To get the continuous case, we replace these with 9c = q, and 9 ^ 2k ~ q, so we have the distribution 

1 (r + /-l)! 



f (r-l)!(;-l)! 9^c+'-' 
{r + l- 1)1 i9c + 9y-^{9c~9y-^ 



{r- 1)1(1 -1)1 (20e)'-+'-i 

Then, if we have / flips, the distribution is 

hf+i^f+l{9) = jTET^ ''%ll]'^ for odd / 

df-i{9) for even /. 

■ 

Lemma II. 2 We let E[g(9)]f represent the expected value of an function g{9) given that f flips occurred. If g{9) is 
odd, then E[g{9)] = 0. If g{9) is even, then 



Elgmi = { 



gidc) forf^O 

-j^fl-—^gf{9c) for odd f ^ (3) 
E[g{e)]f-i for even / > 



where 



9(0) 



e 



9-m = "-Y g^+2{6)= ] Hn{0)d9. (4) 

It then follows from Eq. [I] that the total expected value of an even function g{9) can be written as a sum over the 
E[gi9)]f for odd f as 

\2n4-l \2n+2 \ 

..(^.) + E/[^W]-i((^ + (^)j (5) 
Proof: If / is odd, then by integration by parts, the expected value of g{9) is 

Ebmf = r df{9)g{9)d9 = T {9l - 9^)"^ 9g^,{9)de 

= (/ - 1) ... r (^c - e')'^9g^{9)d9 



(20^)/(/_i!)2 

= (/ - - 3) . . . 2 ^^^^-i / " 0gf^2(9)d9 

^ ' {29c)H^\)^ J-e^ ^ 

/! g/(gc)-ff/Mc) _ , 9f{0c)-gf{-ec) 

2^9ii^\ 2 -^^^ '^••■^ 29l 
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where g-i{d) — and (7,1+2(6') = f 9gn{6)d9. This is zero for odd functions, and we obtain the desired resuh for 
even functions. ■ 

Applying this to 17(6') = a;™ for even m and odd /, we get 

lx3---/ ^„ lx3---(m-l) 



\f 



■ (m + l)(m + 3)---(TO + /) (/ + 2)(/ + 4) •••(/ + /«)■ 

For m — 2, ^[^^]/ — J^- Applying this to Eq. [5l the variance is 

a' = E[e'] = el{\ + '^^) = Ahr. + ^(e't - 1) (6) 

The variance is additive with independent noise sources. For t << Tc, a 
and the A^tr^ part is dominant. 

Theorem II. 3 i?[sin(ma;)] — 0, and 



E[co8imx)] = (i + ^)e'^-'^^^ + - (7) 



where 



v = ^l- (mAT,)2. 

Proof: For the Fourier transform e*™^, we get 

gf{9) = m~^0P^jf_i{m0), (8) 
where jn{x) are the Spherical Bessel Functions of the first kind [1], which are written as the series 

s— 

Then, the expected value for odd / or even / + 1 flips is 

where z = mOc = mtA. From Eq. [51 it follows that 

E\e-^] e-(cos(z) + V i^ttDim^J^ + 



e-(cos(.) + + ^(S)^- 

fc=l ■ fc=0 



Now since zj^i{z) — cos(z), and jn{z) satisfies the equations 

- COS - 2zt = V"— j„_i(z) 



n=0 



we have 



incVz^-2zt = -rin(^), 

— ' n! 

n=0 

£:[cos(m6')] = £;[e""^] = e-^(cos - + Asinc^z^ - A^) 
= e"^ f cos MmtAf ~ {- f + —sine J (mt Ay - (-)2 



= e |^cos(iy TO-^a^ H sinc{t ^ m'^ A'^ 

I cosh(— + sinhf— w) ) = e^^ (cosh(Aw) + v^^ sinhfAw)) 

\ Tc Tc J ' 
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a. Different limits For small Tc we have Atc 1, and therefore 

E[co8{me)]^ e-i'""'^'^". (11) 

If uiAtc ^ 1, then E[cos{md)] « e~'''cos(u). If in addition, Tc ^ 1, then E[cos{m9)] « cos{mdc)- 
Up to first order in A, we have 

i?[cos(m0)] « cos z + A(sinc z — cos z) « cos z + A — . (12) 

3 



A. Multiple sources and Gaussian noise 

If we have multiple independent sources of Random Telegraph Noise, then since e™'^ is a Characteristic Function, 
the expectation value is the product of the expectation value for each source. 



i;[e™^] :=.e-^^^'Y[icos^imtxo^)^-X^ + X^smc^imtxo,)'^-X^] (13) 

i ^ ^ 

If the number of flips A is large, then by the Central limit theorem, the distribution of Q will approach the Gaussian 
(also called Normal) distribution 

P[x)^^j=e-^ . (14) 
(TV27r 

For a single source of noise, from Eq. [51 for large A, cr^ « ^ = b^T^t. This is also (Atc)^A, as one would expect 
from a random walk. 

If we have r telegraph noise sources, with a mean correlation time of r„ and a distribution of telegraph strength 
with mean Am and standard deviation A^., then a « ^■(A^j + A^)Tmi. 

h. 1/ f noise Suppose we have r 1// telegraph noise sources with a frequency in the range [/a, fb\- The distribu- 
tion for / is /(iog(/,)-iog(/„)) 3- ^i'^ce = j, then A = ft, so we have a distribution in A of x{\os(\i)-\og,{\^)) ^ where 
Ai ~ fit. From Eq. [6l it follows that 

/•Ai, 2^ e^^^ — 1 1 
"' = '1 ''^A + -^A^)A(log(A,)-log(A.))''^ 
r6'2 /-Ab -2A _ 2 



log(A6)-log(A,)7^^ ^A2 2A3 
rel fl 1 , -^^(2A-1)_^^(_2A) 



I Ab 



log(Af,)-log(Aa) V4A2 A 4A2 
where Ei [x) is the exponential integral. Now if the A^ are large, the j part dominates, and so 



J2 



a 



A(log(Ab) - log(Aa)) log(r„) - log(rfc) log(^) ' 



= rtA2-^, (15) 



where = j-. This could also be found directly from cr^ « ^ for a single source. If w = ^ = then ct^ — 
rtA?Tb la^l,) ■ Since Tq > t^, HA'^ti, < cr^ < rtA^Ta- If r = 1, and Tq = rt = Tc, this gives the single source result. 

Sometimes by 1// noise, we mean that the noise has a power spectrum of for some a > 1. In this case, then 
by a similar calculation to that of [2], on the domain [0,t;,], a power spectrum of -j^ would give a density of random 

telegraph correlation times of g{Tc)dTc = — drc, and so = E[tc\ — JJ' Tc — — ^^Tb. 

'''h '''b " 

c. Expected values Since P{x) is even, ~ for odd n. For even n, we have 

1 /"OO ^2 

E[x-] = —= / x-e-^ 

1 2 /"OO 2 

' ,2„n-l„--^ loo I / ^\„n-2 2, 



{-a'x-''e'^r^+ (n-\)x--'a'e'^) 
= (n - l)cr2£;[2;"-2] ^ ^"(1 X 3 • • • (n - 1)) : 



2^|! 
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Now £'[sin(TOa;)] = 0, and 



(2n)! ^ (2n)!2"n! ^ n! 

In the case where we have just one flip time Tc, we get the same result as in Eq. 1111 

III. CONTROL 

In this section, we assume that the noise is only generated from some Hamiltonian that can be switched on and off. 
We assume that we need to apply this Hamiltonian for time t. 

Lemma [11.31 gives the expected values of the Fourier functions without control, assuming that the telegraph starts 
off in either state with equal probability. Up to second order in t, this is 

E\cosm0] w 1 w 1 - .1^ 

^ ^ 2 2 

Now, we have the following two possible methods to reduce errors. 

A. Waiting method 

Suppose we wait for a time much greater than the correlation time t' >> r^. This will randomize which direction 
the telegraph is going before we apply a Hamiltonian. If we do this n times, then 

E[cos{me)] = (^(1 + v-^)e^^-^+''^ + (1 - v-^)e^'^-^-^^)'' 

mA— , 



2 ' 2n ' 

and E[sin{Tn9)] = 0, so by applying the Hamiltonian for time n different times, and waiting a while in between, the 
rate of error p is changed to ^ . 

B. Suppressing errors method 

Suppose that we have a Hamiltonian H that we use to create the quantum gate e**^, and it generates a noise 
Hamiltonian TV, so that errors are of the form e*^^ where 9 follows a random telegraph noise distribution. If H and 
N commute, and we can apply a gate R with a low rate of errors that commutes with H and anti-commutes with iV, 
then instead of applying e**'^, we apply the gate 

This breaks the time t into n intervals of equal length, with the direction of the telegraph reversed in between. 
Proposition III.l The expected value of E[sin{m9)] = 0, and 

[cos^mw)] w 1 



n^Tc 

2A„ 



In this case a rate of error of p is changed to j^P- 
Proof: Assume that the telegraph starts going off to the right. Then from Eq. [22l since c — im, 

Ele'""^] = e-^(cosh(A?;) + {- + ^) sinh(Av), 

V AV 



where v = -^/l — (mArc)^. Up to third order in t, this is 



E[e-\ « (1 - A + Ia^ - + ^) + + ^)((A.) + ^)). 
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Up to second order, this is 

E[e"^% « 1 - ^—^X^ + ^mt^{l - A) = 1 - Il^^M! + ^mt^(l - -) 
2 2 Tc 

Now suppose that we have two intervals of time t of telegraph noise, with the direction of the telegraph noise reversed 
in between. The expectation value of e""^ for the first is E[e""'^]r. Reversing the initial direction sends 9 —0 for 
the second, giving £'[e~""^]r, and so the total expected value of e™^ is 

S[e™''] = E[e'-^%E[e~'-^% « (1 - l!!^M!)2 + (^iA(l - - f 

2 Tc 

w 1 - 27n^l^^T-^f 

For ^ pairs of intervals with ^ time per interval, this becomes 



From Eq. [22l if the telegraph starts off to the right, instead of i?[sin(mx)] = 0, we now have 

E[sm[mx)\ — e sinhlitj 



mA 



Up to first order in t, this is mAt(l — 



imh{t Vtc'^ - m^A^) . 



C . D r awbacks 



Both methods assume that the telegraph is equally likely to go in either direction. The waiting method assumes 
that general decoherence isn't a problem. The suppression method assumes that the gate R can be quickly performed 
without errors. These assumptions are unrealistic, and other types of noise will likely be created, but these methods 
could allow for a significant reduction of the magnitude of the noise. 



IV. EXAMPLE 



From a superconducting qubit system ^ we have the noise 



je(z®i+i®z) 



where 9 has a random telegraph distribution. Now, 



E ^ cos^ 91 ® I + isa\9 cos9{I Z + Z I) - sin^ 9Z (g> Z, 



(16) 



so the probability of no Z error no = cos"* 0, each of the probabilities of a Z error on exactly one qubit are n\ — 
sin^ 9 cos^ 6*, and the probability of a Z error on both qubits is — sin* 9. We have that 



no 
ni 
n2 



E[cofi* 9] 
E[svc? 9 cos^ I 
E[s\n* 9] 



"3 1 

8 2 


1 " 

\ 




E[l] 


1 


1^ 




E[cos29] 


3 1 
-8 2 


8 . 




£; [cos 461] 



where the _B[cos(mx)] are given in Thm.. III. 31 Then from equation II 2[ we have that up to 1st order in A and 3nd 
order in 9c, 



no 

ni 
n2 



cos 



sin2(e'2)cos2(6'2) 



sm 



-X9i 



7 





3 


1 


-2a^ 


no a 


' 8 " 




ni ft 


1 


1 




' 8 " 






3 


1 




n2 a 






' 8 " 


-r 





d. Gaussian noise Suppose we assume the noise is Gaussian distributed with a standard deviation of cr, then 

1 3 11^ , . 

r -8-2^+r' 

where r = e~^"'^ . This gives, up to cr^, np w 1 — 2o-^ + 5a^, ni fa — 4cr^, n2 ~ 3cr^. 

e. Control The system Q has a 2 qubit Hamihonian H = X iSi X + Y Y , which commutes with the error 
Hamiltonian oiEH = I®Z + Z®I. The gate 

R^X®X 

commutes with iJ, and anti-commutes with Eh- This is very useful, because R is entirely composed of local gates, 
which have a much lower rate of errors. This is similar to quantum "bang-bang" control 

V. GENERAL DERIVATION 

In section ini we considered the expectations of a 2 state random telegraph source with a correlation time tc- In this 
section we consider that the correlation time depends on which state we're in, that is we have chance of flipping 

per unit time if the telegraph is in the positive state, and ^ of flipping time if the telegraph is in the negative state. 
This is relevant to a physical system at finite temperatures. 

Also we assume that the telegraph starts off in the positive state. Note that if for the positive state i?[e™^]+ ~ 
w(m. To, Ti), and there is probability of pi of starting in the positive state, and po of starting in the negative state, 
then 

£;[e^"«)] =Po£^[e""(-«)]_+Pii?[e™%=po«(-™,ri,To)+piz;(TO,To,ri). 

Lemma V.l Suppose the telegraph starts off in the positive state (so that the noise parameter is increasing). The 
distribution of is given by a sum over all of the possible number of flips 

oo 
/=0 



where 



h{9)S{e-e,) tff^o 



df{0) = lmj^^}^Plf^{Oc + OM-e^)i-^ for even f>0 (20) 

[mfi^^^r^iO^. 0^)^ for odd f 

h{e) = We-^i w (21) 

on the domain [— 0c,^c]j where 6c — tA, A,; = 

Proof: There are no flips with probability e^^^, in which case 6 = 6c, which gives do{9) = e~^^5{6 — 9c). 

Suppose we have r intervals where the telegraph is in the positive state and so the noise parameter is increasing, and 
I intervals where the telegraph noise is in the negative state and so the noise parameter is decreasing, and we have q 
steps. Then we have g{r, k)g{l, q—k) ways to end up at position 2k — q. Suppose we start off in the increasing telegraph 
state, and the probability of flipping if we're increasing is t^^ , and if we're decreasing is Tq^. The probability of being 
at position j afterwards is the coefficient of in 

V-(l - r^'f-'{r^'y{l ror-'-^+\r,'r-'g{r, k)g{l,q- k)x''^-'^ 



fc=o 



_<ti^, _i,r_ifc'' ^{q-ky ^ 2k-o 
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We make this continuous, so we have 9 = 2k — q, 6c = q, and pick up a factor i from and so have 



— e e Tq 



(r-l)!(;-l)! 

h(e) (^c + grH^c^g)'-^ 
^ 2^-+i-iT[T^-\r - iy.{i - ly. 



If we have an even number of telegraph flips / > 0, then r=^ + l,Z = ^. If we have an odd number of telegraph 
flips, then r = 1 ~ Plugging these into the previous equation produces the desired result. ■ 

Lemma V.2 

„cenln2 /l2\n I/) _ 2"n!(— 1)" „n \ „ce 



where the pn{x) are the Carlitz Bessel polynomials described in 0/. 
Proof: If we let 

9o{e) = e^' 

9n+l{0) = J 9n 

then by integration by parts. 



2n / Qk+i - 2ne{el - O'^T-^dO - 2"n!(5„+fc+i(0,) - g„+fc+i(-0,)), 



so 

a„ = 2"n!(5„+i(6'c) - 5„+i(-6'c)). 
Now, from the the formula for the Carlitz Bessel functions, it can be shown that they satisfy the differential equation 

P'n+l{x) -_p„+l(x)(l +XPn{x) = 0. 

Note that this implies that Pq{x) = 1. If we let 

p.n+l{-c6) {-IT B 

-e ' 

then this satisfles the differential equations 5^4+1(6*) = 6gn{9), gQ{9) = e'^^, and so we get the desired result. 

■ 

Theorem V.3 The expectation of a characteristic function e"""^ is 

£;[e™«] = e"^^(coshu+ ^ii^sinhw) (22) 

u 

where 



^c^9l + AoAi - tJc^^^ + V Vf 1, (23) 



w/iere c = + im. 



Proof: From Lemma IV.ll the contribution from an odd number of flips / is 



71=0 

By integration by parts, 

2n 



c 



and so from Lemma IV.21 if z = c6c 



£;[e™«]odd = -Ai^— ^5](-^)"-(p„+i(-z)e^+p„+i(z)e-^) 



Now, by differentiating the formula 



from 3, we get 



so if t = 



kl 

k=0 



;(l-\/l-2t) 



aq+ai sinh(zA/l — 2t') , aq+ai sinhu 
= Aie 2 , — - = Aie 2 



For even / > 0, the contribution is 



and 



so if z = C^r 



In the sum, if n = 0, 



■^0 ~^ 1 00 

?l— 1 



{{pQ{-z)e +po{z)e ) 



Aq + Ai 
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so we can include / = 0, and so we get a sum over all even /, with t = — ^p^, 



^[e Jc 



e 2 ^-^ t 



2 ^ n 

Tl = 



Aq+Ai oo , 

e 2 . ^ t 



2 ^ n 

2 u 

= e 2 (cosnuH — smliTi) 
u 
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